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Abstract. We expand a discrete-time lattice sine-Gordon equation on mul- 
tiple lattices and obtain the partial difference equation which governs its far 
' field behaviour. Such reduction allow us to obtain a new completely discrete 

CNJ I nonlinear Schroedinger (NLS) type equation. 
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1. Introduction 



Reductive perturbation technique [IHIIII] has proved to be an important tool to 
^ . find approximate solutions for many important physical problems by reducing the 

given nonlinear partial differential equation in the far field often to an integrable 
^ I one. Recently this approach has been extended to the case of equations living on 

lattices [Ul |H] . Here we apply it to the case of a discrete-time lattice sine-Gordon 
^ ■ equation. 
O ! In Section |21 we briefly describe the discrete perturbation technique and in 

Section El we apply it to a discrete-time lattice sine-Gordon equation. 



2. The discrete perturbation technique 



■ The aim of this section is to flx the notation and to introduce the formulae 

^ ! necessary to reduce lattice equations in the framework of the discrete reductive 

perturbation technique [HI IHI • 

Given a lattice deflned by a constant spacing h, we will denote by n the running 
index of the points separated by h. In correspondence with the lattice variable 
5^ ! n, we can introduce the real variables x = hn. 

We can deflne on the same lattice a set of slow varying variables by introducing 
a large integer number N, deflning a small parameter e = A^~^ and requiring that 

(1) Uj = n. 

This correspond to sampling points from the original lattice which are situated 
at a distance of h between them. If we set them on a lattice of spacing h, the 
corresponding slowly varying real variables Xj are related to the variable x by the 
equation Xj = x. 

Let us consider a function f = fn deflned on the points of a lattice variable 
n and let us assume that /„ = gni,n2,...,nKi / depends on a flnite number K 
of slow varying lattice variables j = 1, 2, . . . , deflned as in eq. ((T)). We are 
looking for explicit expressions for, say, /„+i in terms of gni,n2,...,nK evaluated on 
the points of the ni, n2, . . . , rix lattices. At flrst let us consider the case, studied 
by Jordan jl], when we have only two different lattices, i.e. K = 1. Using 
the results obtained in this case we will then consider the case corresponding to 
K = 2. The general case will than be obvious. 
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1) K = 1 {fn = fi'm)- In Jordan j;4 we find the following formula: 

JL /k\ ^ /rl 

(2) A'^g^,=J2i-^)''H . Un.+. = 5^-7P(^,A:)A7n. 

i=0 ^ ^ i=k 

Here the coefficients P{i, k) are given by 



(3) P{^,k) = J2^^Sr& 



k 

k 



where uj is the ratio of the increment in the lattice variable n with respect to that 
of variable rii. In this case, taking into account eq. u = N. The coefficients 

and are the Stirling numbers of the first and second kind respectively jSj- 
Formula Q allow us to express a difference of order k in the lattice variable ni 
in terms of an infinite number of differences on the lattice variable n. The result 

can be inverted and we get: 

°° I.} 

(4) A^/„ = ^-Q(^,fc)A^^„„ 

i=k 

where the coefficients Q{i, k) are given by with u = = e. 

To get from eq. © a finite approximation of the variation of = /„ we need 
to truncate the expansion in the r.h.s. by requiring a slow varying condition for 
the function /„. Let us introduce the following definition: 

Definition. /„ is a slow varying function of order p iff A^"*"^ /n = 0. 

From the above definition it follows that a slow varying function of order p is 
a polynomial of degree p in n. From eq. © we see that if is a slow varying 
function of order p then Ap~^^ = 0, namely gm is also of order p. Eq. (jU 
provide us with the formulae for fn+i in terms of gn^^ and its neighboring points 
in the case of slow varying functions of any order p. Let us write down explicitly 
these expressions in the case of g^ of order 1 and 2. 
• p = 1. Formula (jl)) reduces to 



A/„ = J^^9ni; 



i.e. fn+i reads 



fn+l = gm + -gm)+0{N 2) 



p = 2. From eq. (0} we get 



1a 1-^a2 

A/„ = ^A^„^ + ^^A ^„^, 
and thus fn+i reads 

fn+i = gm + ^(-fi'„i+2 + 45-^1+1 - 3 5f„J + 
(5) +^{gn,+2-2gn,+i+gn,)+OiN-'') 
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In the next section we will consider the reduction of an integrable lattice sine- 
Gordon equation. It is known |12j that a scalar differential-difference equation 
can possess higher conservation laws and thus be integrable only if it depends 
symmetrically on the discrete variable. The results contained in (0)) do not provide 
us with symmetric formulae. To get symmetric formulae we take into account 
the following remarks: 

(1) Formula Q holds also if h is negative; 

(2) For a slow varying function of order p, we have /SP = /SP fn+i, for all 

When fn is a slow varying function of odd order we are not able to construct 
completely symmetric derivatives using just an odd number of points centered 
around the ni point and thus fn±i can never be expressed in a symmetric form. 

From eq. using the above remarks we can easily construct the symmetric 
version of (0). We get: 

(6) fn+i = + ^(^m+i - gn,~i) + ^^(^m+i - 2^ni + g^-i) + 0{N'^). 



II) K = 2 {fn = gnx,n2)- The derivation of the formulae in this case is done in 
the same spirit as for the symmetric expansion presented above, resulting in eq. 
(P). Let us just consider the case when p = 2, as this is the lowest value of p for 
which we can consider /„ as a function of the two scales rii and n2. From eq. (jH) 
we get: 

(7) gn,+l,n, = gn,,n, + NA,U,n + ^NiN - l)AlU,n, 

(8) gnun2+l = gm,n, + A^fn.n + ]^N\N^ - 1)A2/„,„. 

Here the symbols Ai and A2 denote difference operators which acts on the first 
and respectively on the second index of the function /„^„ = gni,n2, e.g. Ai/„^„ = 

fn+l,n fn,n- and A2fn,n fn,n+l fn,n- 

Let us now consider a function gni,n2 where one shifts both indices by 1. From 
eq. (0), taking into account that, from eq. (H)), for example, gni+i,n2 = fn+N,n, 
one has: 

(9) S'ni+l, 712+1 5'rai,n2+l 

Using the result (jH)) we can write eq. © as 

+iVAi/„,„ + iV^AiAa/n,. + ^iV3(iV2 _ l)A^Alf^^n + 
+liV(iV - l)A?/„,„ + N\N - l)A?A2/„,„ + 

(10) +\n\N' -1){N -I) AlAlU,n. 

As, using the second remark, the second difference of depends just on its 
nearest neighboring points, the right hand side of eq. depends, apart from 

fn,n gn\,n2) fn,n+lj fn,n—l} fn+l,ni fn—l,ni fn+l,n+li fn+l,n—l^ fn—l,n+lj and 

fn-i,n-i, i-e. 8 unknowns. Starting from eqs. ((Tj), (jSj) and (jTU)) we can write 
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down 8 equations, using the first remark, which define gni+i,n2} 9ni-i,n2j 9ni,n2+ii 
9nun2-i, 9ni+i,n2+i, fl'ni+i.na-i, 9ni-i,n2+i, and gni-i,n2-i temis of the fuuctious 
fn+i,n+j with = 0, ±1. Inverting this system of equations we get fn±i in 

term of 9ni,n2 aiid its shifted values: 

/nil ~ 9ni,n2 i 2 A'^ i9ni+l,n2 fl'ni — 1,712) 2 ]\f2 i9ni,n2+l 9ni,n2~l) ~^ 
1 

{9ni+l,n2 9n\,n2 ~^ 9ni — l,n2) ~^ 



2 m 



+ ^ ^3 (5'ni+l,n2+l ~ 5'ni-l,n2+l ~ 5'ni+l,n2-l + 9ni-l,n2-l) + 

(11) +0{N-^). 

It is worthwhile to notice that the two lowest order (in A^~^) terms of the expan- 
sion pi|l are just the sum of the first symmetric differences of Qm and gn2- Thus 
in the continuous limit, when we divide by h and send h to zero in such a way 
that X = hn, xi = hrii and X2 = hn2 be finite, we will have f^x = ^9xi + ^'^9x2- 
Extra terms appear at the order A^~^ and contain shifts in both rii and n2. 

When fn is a slow varying function of order 2 in rii it can also be of order 1 in 
712. In such a case eq. (jHl) is given by 

(12) 9ni,n2+l = 9ni,n2 + ^^2/n,n- 

Starting from eqs. ((Tj), (fT^ and a modified (fTUj) we can get a set of 8 equations 
which allows us to get fn±i in terms of 9ni,n2 its shifted values. In such a 
case fn±i reads 

/nil 9ni,n2 ^ 2 A'^ i9ni+l,n2 9ni-l,n2) ~^ 'J^{9ni,n2±l 9ni,n2) ~^ 

I 

(13) +^-^(^ni+l,n2 - 2 9m,n2 + ^ni-l,n2) + 0(A^"^). 

It is possible to introduce two parameters in the definition of rii, n2 in terms of 
n. Let us define rii = {nMi)/N,n2 = {nM2)/N'^, where Mi and M2 are divisors 
of and A^^ so that rii and n2 are integers numbers. In such a case eq. pi|) 
reads 

f - -^Mif _ \ + ^( _ 

/nil 9ni,n2 ^ + 9ni — l,n2) ^ 2 ^^"1 '"2+1 9ni,n2~l) ' 

Ml ^ 

2 \9ni+^,n2 9ni,n2 ~^ 5'rii — 1,712) 

M1M2 

+ ^jyS (fi'ni+1,712+1 ~ 5'ni-l,n2+l ~ 5'ni+l,n2-l + 5'ni-l,7i2-l) + 

+0{N-^) 

and eq. (fT^ accordingly. 

Let us consider the case of two independent lattices and a function ^ defined 
on them. As the two lattices are independent the formulae presented above apply 
independently on each of the lattice variables. So, for instance, the variation 
fn+i,m when the function fn,m is a slowly varying function of order 2 of a lattice 



MULTISCALE EXPANSION OF THE LATTICE SINE-GORDON EQUATION 5 

variable rii reads 

fn+l,m 9n\,m ~\~ 2 jV ^3n-i+l,m 9ni — l,m) ~l~ 

+ ^^(5'ni+l,m - 2 gn^^m + 9n^-l,m) + 0(iV"^). 

A slightly less obvious situation appears when we consider fn+i,m+i^ as new 
terms will appear, see Levi et al. for the formulae in this case. 

3. Reduction of the lattice sine-Gordon equation 

A discrete analogue of the sine-Gordon equation is given by the following 
nonlinear PAE ^ 

1 4 

(14) Mn+l,m+l = _ 4 • 

This equation involves just four points which lay on two orthogonal infinite lattices 
and are the vertices of an elementary square. When written in polynomial form 
has quartic nonlinearity. In eq. ()14j) M„,rra is the dynamical (real) field variable at 
site (m, n) G Z X Z and p,q E M. are the lattice parameters. These are assumed 
different from zero and will go to zero in the continuous limit so as to get the 
continuous sine-Gordon equation. 

To get a nonlinear dispersion relation we carry out the change of variable 
Un,rn ^ p/q + Un,m- The linear part of the resulting equation is given by: 

(15) icr-l){u ) + (a + l)(M„+l,m + Mn,m+l) = 0, 

where a = p^q^. 

The general solution of eq. (jl5j) is written as a superposition of linear waves 

En,m = exp[i(/cn — C(j(A;)m)] = z"'Q"^. The dispersion relation for these linear 
waves is given by 



(16) n = e-= (- + 1)^ + ^-1 



lui 



namely 



iv = — arctan 



a-l)z + a + l' 
2 a sin k 



(cr2 + 1) cos A; + _ 1 
From eq. ()l(ij] . by differentiation with respect to k, we get the group velocity: 

Aa z 

" ~[{a+l)z + a~l][{a-l)z + a + iy 

We now look for real solutions of the nonlinear equation ()14j) written as a 
combination of modulated waves: 

oo oo 

(18) Un,m = J2 ^""^^nl iEn,mr + Yl ^""^^^1 (^".-)^ 

s=0 s=l 
(s) (s) 

where the functions ipn/m are slowly varying functions on the lattice, i.e. ipn/m = 
4'n},m-i,m2 and = N'-^. By b we mean the complex conjugate of a complex 
quantity b so that, for example, En^m = {En^m)~^- The positive numbers Ps are 
to be determined in such a way that : 

(1) /?i < /?s V s = 0, 2, 3, . . . , oo. In general it is possible to set Pi = 1. 
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(2) In the equation for ijjnjn = ipn,m, the lowest order nonhnear terms should 
match the slow time derivative of the linear part after having solved all 
linear equations. This will provide a relation between 7 and the f^g. 

Introducing the expansion (fTHjl in the PAE obtained from equation (fT^ after 
the change of variable Un,m ^ p/(l + Un,m, we analize the coefficients of the various 
harmonics [En^mY for s = 1, s = 2 and s = and, as assuming that (3s increases 
with s, the nonlinear terms will depend only on the lowest s terms, we came to 
the conclusion that we can choose 7 = 1, /3o = 2, /S^ = s, s > 1. The discrete slow 
varying variables ni, mi and 1712 are defined in terms of n and m by: 

Ml n M2 m n 

For s = 1 we get, at the lowest order in e, a linear equation which is identically 
solved by the dispersion relation (fTT)|l . 
At we get the linear equation 

(19) + M2ll [(a - 1) 2 + Cr + 1] (^„,,m,+l,m2 - V'ni.mi-l.ma) = 0, 

whose solution is given by 4'ni,mi,m2 = 0n2,m2? n2 = rii — rrii, provided that the 
integers Mi and M2 are choosen as 

(20) Mi = Sn[{a-l)z + a + l], M2 = S z [{a ~ 1) n + a + 1], 

where 5* G C is a constant. Defining 5* = pexp(i6'), p G IR+ and — vr < < tt, we 
can choose 6 and p in such a way that Mi is an integer number: 

(cr + 1) sin k 



9 = — arctan 



(cr + 1) cos k + (7 — 1 



+ i7l, £GZ, 



p = (-l)'Mi — [{a^ - 1) cosA; + + 1]-^/^. 
V2 

The request that also M2 is an integer impose a constraint on k as M2 = cj^^Mi, 
i.e. u^k G Q, see eqs. (ITTjl and (j^ . Let us notice that also n2 = ni + mi solves 
eq. (|Tn|l by an appropriate choice of Mi and M2. 

At e'^ we get a nonlinear equation for which depends on '4)n}m2 and 

tyn2,m2 ■ 

4>n2,m2+l ~ 4>n2,m2 + Cl (0n2+2,m2 + <Pn2-2,m2 ~ (p '112,1712) + 
~l~ C2 (0^2+1, »n2 ~l~ 0n2 — l,m,2 ^ 0n2,m2) ~^ '-'3 0n2,m2 I0n2,m2 I ~l~ 
(21) + C4 ^i°!m2 '^n2,m2 + C5 ^^^.^^ (pn2,m2 = 0, 

where the Cj's, 1 < < 5, are known coefficients depending on z, S and the lattice 
parameters p, q. 

The functions ipn2,m2 and tpn2,m2 that appear in equation (j7T|) are obtained 
considering the equations for the harmonics s = 0, at the third order in e, and 
s = 2 at the second one. We get: 



r22) =-\d |2 ib^^^ =--(d ? 

V^^/ Vn2,m2 \Yn2,m2\ ? rra2,m2 o lrr!,2,m2j • 

P P 
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Inserting ipn2,m2 and il)n2,m2 given by eqs. in equation (|7T|) we obtain the 

following nonlinear lattice equation: 

i (0n2,m2+l ~ 0n2,m2) + Cl {4>n2+2,m2 + '^n2-2,m2 ~ 2 4>n2,m2) + 
(23) + C2 (0„2+l,m2 + 0n2-l,m2 ~ 2</)n2,m2) + C3 0n2,m2 I0n2,m2 P = 0; 

where the coefficients q's, 1 < < 3, can be computed also recalling that z = 
exp(ifc). They read: 

. M|(o-- 1)[((T + l)(cosA; + isinA;) + cr + 1] 



16 a 

M| ((T - 1) [(a + 1) cos /c + (T + 1] 

la ' 
2 (a^ - 1) sin^ A; 



Cl = 

(24) C2 = 

" [((T2-l)cosA; + a2 + l)]' 

The coefficients (E^ depend just on the integer constant M2. 

Taking into account that 0n2,m2 is a slow varying function of order 2 in n2 we 
can, using the remark |21 substitute 4>n2+2,m2 by 3 0„2+i,m2 ~ 3(/)„2,m2 + 0n2-i,m2 
and 0„2-2,m2 by 3 0n2-i,m2 - 3 0„2,m2 + 4>n2+i,m2- In such a way eq. (j2SI) becomes 

i (0n2,m2+l ~ 0n2,m2) + (4Cl + C2) {4>n2+l,m2 + 0n2-l,m2 ~ 2 0„2,m2) + 

(25) + C3 0^2,^2 |0„2,m2 1^ = 0, 

where 

_ „ M|((t2 - 1) sin A; 

(26) 4ci + C2 = ^ . 

4 cr 

The PAE is a completely discrete and local NLS equation depending just 
on neighboring lattice points. At difference from the Ablowitz and Ladik pP 
discrete NLS, the nonlinear term in equation ()25|1 is completely local. The PAE 
()25|1 has a natural continuous limit when 1712 — > 00 and 722 ^ 00 which, as the 
coefficients (|26p and C3 are real, is just the well known integrable NLS equation. 
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